In-situ monitoring and control of mechanical wavefunction collapse: an 

optomechanical scheme 



A. A. Ganga|] 

ARC Centre for Engineered Quantum Systems, School of Mathematics and Physics, 
The University of Queensland, St Lucia, QLD 4 072, Australia 
(Dated: XX March 2013) 

The optomechanical system demonstrated by Sankey et al. [Nat. Phys. 6, 707 (2010)] enables 
different types of parametric couplings between optical fields and mechanical motion to be achieved 
simultaneously, including one that is purely quartic to lowest order in the mechanical motion. We 
propose and analyze a scheme that uses this capability to monitor and control in-situ the conditional 
phonon number (nb) variance of a mechanical mode via continuous weak measurements of nb and n\ 
while simultaneously cooling to {nb) ~ 1. Simulations show a monotonic decrease in the steady-state 
phonon number variance with increasing measurement strength. 
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Introduction. — Creating, manipulating, and observing 
quantum phenomena in macroscopic mechanical degrees 
of freedom is of both fundamental and applied inter- 
est, and a large amount of activity toward these ends 
has taken place within the context of optomechanics pQ , 
wherein optical field modes are coupled to the motion of 
mechanical resonators. Mechanical energy quantization 
[2H1], mechanical mode entanglement [5], and detection 
of space-time quantization effects [5J, for example, have 
all been considered in such a setting. Particularly well- 
suited to mechanical energy quantization detection is the 
" membrane- in-the- middle" optomechanical system [7] , in 
which a dielectric membrane is suspended in the middle 
of an optical cavity and orthogonal to the cavity axis. 
Depending on the equilibrium position Xq of the mem- 
brane along the cavity axis, the system can exhibit a 
modulation of the energy of a full cavity optical mode a 
that could be either linear or quadratic in the mechanical 
membrane displacement x along the cavity axis from x$: 
H cx xo) a or H oc x 2 o)a. With the mechanics in the 
quantum regime and under the rotating wave approxi- 
mation, the latter becomes for a single mechanical mode 
b: H oc b^ba^a, thereby providing a channel for continu- 
ous QND measurement of mechanical energy by driving a 
and continuously monitoring the transmitted signal. The 
linear in x coupling provides a channel for actively cool- 
ing the mechanical mode to low occupation numbers [S] ; 
driving a at the red mechanical sideband induces a net 
up-conversion of the pumped photons via absorption of 
mechanical quanta. Further, tilting the membrane with 
respect to the cavity axis can change the optomechan- 
ical modulation of select optical spectra from x 2 to a; 4 
at lowest order in x [9], providing a channel for QND 
measurements of tfb + (b^b) 2 . Also, an examination of 
the full optical spectra of the system reveals that the x, 
x 2 , and x couplings may all be achieved simultaneously 
with independent optical channels. 

In this paper we show that the capability of simul- 
taneous x, x 2 , and x 4 optomechanical coupling can be 



used to monitor and control in-situ the variance of the 
phonon number distribution of a macroscopic mechani- 
cal mode that is coupled to a thermal bath of very high 
average occupation number while simultaneously actively 
cooling the mechanical mode to the single quanta regime. 
The essential physics behind this scheme is the interplay 
between narrowing of the number distribution due to ac- 
quisition of information via measurement and broaden- 
ing of the distribution due to loss of information to the 
unmonitored dissipation channels. The continuous weak 
measurement of n& through the x 2 coupling effects the ac- 
tion of collapsing the phonon distribution toward a single 
(random) Fock state, while coupling to dissipative chan- 
nels induces a broadening of the distribution toward a 
thermal state. The (time- averaged) steady-state between 
these competing processes has a finite phonon number 
variance. Increasing the n\> measurement strength, which 
may be done in-situ by increasing the drive strength 
on the optical mode coupled to x 2 , results in a smaller 
steady-state phonon number variance because the infor- 
mation acquisition rate is thereby increased. Simultane- 
ously, the information from the x A measurement channel 
may be combined with that from the x 2 channel to deter- 
mine the steady-state variance. Thus, the wavefunction 
width of a macroscopic mechanical mode may be moni- 
tored and controllably collapsed in the number state ba- 
sis. The simultaneous mechanical mode cooling through 
the x coupling serves the purpose of lowering the effective 
bath temperature of the mode, thereby facilitating the 
viability of the scheme at smaller x 2 coupling strengths 
than would be required otherwise. 

Model. — We first follow some of the analysis of [H [10] 
for the case of x 2 optical spectrum modulation in the 
case of a two-sided cavity. In the membrane-in-the- 
middle system, when the membrane is orthogonal to 
the cavity axis the finite optical transmittance of the 
membrane and the finite optomechanical coupling give 
rise to the following Hamiltonian at select values of 
Xq and valid for small x of a single mechanical mode: 
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Hi = Tiu)ia\ L aiL + hu)ia\ R ai R - hJi(a[ L a 1R + H.c.) - 
Ttgi(x/x Z p{)(a\ L aiL — a\ R am), where Ji is proportional 
to the transmittance of the membrane for clil/ir, 9i is an 
optomechanical coupling constant, x zp f is the zero point 
motion of the mechanical mode, and clil/ir are one of 
the cavity mode pairs, degenerate in frequency but of 
different spatial mode functions, that would arise on the 
left/right of the membrane with frequency U)i if there 
were no interactions (Ji = gi = 0). In terms of the 
full cavity modes a x _± = (aiL ± aifl)/i/2, the Hamilto- 
nian is Hi = Hi + H[ mt \ where H[ ' = ?iw 1 _a{_a 1 _ + 

hu)i+a\ + ai + , H{ nt > = -hgi(x/x zpi )(a\ + ai_ + a|_a 1+ ), 
and uji± = U)i =p J\. Thus, though ai.± are degenerate 
when the interactions are off, turning on the interactions 
alters the optical spectrum of ai t ± such that, in the phys- 
ically relevant case of Ji ^> gi(x / x zp f) , the transmittance 
lifts the degeneracy by an amount 2hJi, and the optome- 
chanical interaction induces a further perturbation of the 
spectrum that is quadratic in x for small x. 

As the spectral splitting of two modes due to a lin- 
ear coupling between them is always quadratic in the 
coupling strength, we can go beyond [5J [TU] to analo- 
gously model the case of the x 4 spectrum modulation 
of the full cavity optical modes a 2y ± = (a 2 L ± o*2r)/V2 
that was demonstrated in H 2 = H^ + H^"^, 

where H% — huj 2 -a\_a 2 - + ^2+i2+ a 2+) H 2 nt ^ = 
-hg 2 {x/x zp i) 2 {a\ + a 2 - + <4_a 2+ ), and ui 2 ± — ^2 T J%- 

We are interested here in the case of simultaneous 
x, x 2 , and x A coupling for the fundamental mechan- 
ical mode b so that the full Hamiltonian is given by 

H = hfltfb + Hi + H 2 + i3l,drive + #2, drive + #diss, 

where ft is the fundamental mechanical mode frequency, 
drive = h(e*aj + e luJj+t + H.c.) encapsulates coherent 
drives of amplitude tj on optical modes a,j+, and H d ^ ss 
encapsulates all of the intrinsic and induced dissipation 
channels, further discussed below, for the relevant opti- 
cal and mechanical modes, including the mechanical side- 
band cooling bath due to the x coupling [11] . 

We may now make explicit the different order in gj 
perturbations of the spectra due to i/j mt ' by using 
the approach of [8j; the strategy is to make a time 
dependent canonical transformation determined by the 
O(gj) mechanical dynamics to obtain a transformed 
system Hamiltonian for which the lowest order in the 
perturbation are quadratic: e l ( Sl+S2 \h£ltf b + H\ + 
H 2)e -r(s 1+ s 2 ) = mb ] h + H (o) + H (0) + i^tfM)] + 

l 2 [S 2 ,H { t t] ] + i\S x ,H$ at) ] + |[5 2 , H? nt) ] +0(g? ), where 

Sj = X indJ p/h - Pind.jx/h, X = X zp f(tf + b), p = 

ih(tf — b)/x zp f, x zp f = yjh/2m£l, and m is the me- 
chanical mode effective mass. As mentioned in [S], the 
quadratic optomechanical interaction induces, to first or- 
der in gi, the mechanical motion x- m d,i = Fi(t)/m(Q 2 — 
4J 2 ) and p ind ,i = Fi{t)/(n 2 - 4J 2 ), where Fj(t) = 



^L(2tfe mt + H.c.y-H^+^-e""^* + H - c -)- Similarly, 
the quartic optomechanical interaction produces a radi- 
ation pressure force with two frequency components so 
that, to first order in g 2 , x ind . 2 = | m(n 2_^j 2 -o) 2 ) + 

5 ti>(n»-Qj a +n)») and ^,2 = "ii ind , 2 , where F 2 . ± = 
^(a^a_&t e -»(2J±n)t + jj c ) xhus finding the expan- 
sion of the system Hamiltonian in powers of gj , we plug 
into the full Hamiltonian, transform to a picture mov- 
ing at the zeroth order optical and mechanical dynamics, 
make the rotating wave approximation, and drop the first 
order in gj terms of the expansion (as they do not mod- 
ify the optical spectra) to find the effective Hamiltonian 
H cS = H[ 2 cS + H ( 2 2 cS + H' ldlWe + H 2drWc + H diss + 0{g Z j), 
where 

+ h % (uf-Q?) ( 2ni - ni+ + ni+ + ni -) ' 

(1) 

2 

H Ss = h Y ( U2 - ~ n2 +) { A2n b + AlUb + As ) 

- h 92 j2^Q2 n 2+^2- (3 + 4n b ) , (2) 

n b = tfb, H' jATivc = h(e*a j+ + H.c), H diss is unchanged 
as it is modeled with Lindbladian superoperators (see 
below) that are invariant under transformation to the 
zeroth order dynamics, A 2 = j- + jtz^s, Ai = j- + 

a 3 = n(3+n/j 2 )/(j 2 -n 2 ) + a/J2+2J2/(k- 
n 2 )), and a 4 = n(3-n/j 2 )/(j 2 -n 2 )+(i/j 2 +2j 2 /(j 2 - 

fi 2 )). These expressions are valid for all values of the 
ratios Cl/Jj. We see here that the coupling enhancement 
when 2 Ji — ft <C Ji, Q noted in [5j for x 2 coupling has an 

analogous counterpart in the case of x 4 coupling when 

(2) 

J2 — it <?; J2, fi; if Jj 3> " the second lines of H^ o ; ff and 
(2) 

H\ c ; ff become negligible. 

It is clear that the frequency of mode ai,± is sensitive 
to rib and the frequency of a 2 ,± is sensitive to both and 
n 2 . When modes a,j+ are continuously driven, their out- 
puts therefore yield continuous information on nb and n 2 . 
However, in the case that 2Ji — £1 <C Ji,^, the output 
of (Z14. is actually sensitive to (ni_ — rib). Although we 
model the situation where cij- are left undriven, Hj lnt ^ 
support secondary Raman processes whereby quanta 
from a,j + combine with phonons to scatter into a.j-- Let- 
ting Kj± be the total intrinsic dissipation rates of %±, 
reference [S] determines this to occur at a rate 71,^6 = 
= g 2 ni + ni,Ki- / (2Ji — Q) 2 for scattering from 
ai + to ai_, and an analogous golden rule calculation 
yields 72,^ = K 2+yR n 2+ = g 2 n 2+ n b n 2 ^/(2J 2 -2n) 2 for 
scattering from a 2+ to a 2 -. Therefore 71. r is also ampli- 
fied when 2 Ji — f2 <C Ji , il, and a non-negligible ni_ may 
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result. In the case of quantum jump measurements, al- 
though not noted in the study of [8] , this Raman process 
can thereby register as a double jump in the output sig- 
nal of ai+ due to its sensitivity to (n%- — rib). Similarly, 
in the case of J 2 — il <C J 2 , fi, a 2+ becomes sensitive to 
ri2- and J2,R is amplified so that n^- may become non- 
negligible. We are interested here, however, in the case 

of long time-averages of the outputs of dj+ and assume 

(2) 

that even in the cases where the second terms of Hj ' 
become significant, their time-averaged contributions to 
the output signals may be subtracted away by, for ex- 
ample, independently monitoring the output of aj-. For 
simplicity, then, we deal with the following model Hamil- 
tonian: 

ffmodel = - 2 G\n 1+ n b - -G 2 n 2+ (n\ + An b ) 

+ -^1, drive + ^2, drive + ^dissi (3) 

where G 1 = 2gf(^ u + ^j^jj), G 2 = g 2 2 A 2 , A = 
A\/A 2 , and dissipation due to the environment and Ra- 
man processes is contained in -ffdiss- 

Collapse monitoring and control. — The protocol for in- 
situ control and monitoring of phonon number wavefunc- 
tion collapse is as follows. As established in [TJ, the 
steady-state phonon number distribution under contin- 
uous measurement of n b is the result of a competition 
between conditional narrowing of the number distribu- 
tion due to the information acquired in the measure- 
ment record and broadening due to loss of information 
to the unmonitored dissipation channels. Although the 
conditional number distribution itself fluctuates in time 
due to the continuous weak measurement, for the case 
of constant measurement and dissipation rates the time- 
average of the conditional number distribution variance 
(in steady-state) is constant. If the dissipation rates 
are constant, increasing the measurement strength on 
n b results in a smaller time-averaged steady-state condi- 
tional variance for n b . As the n b measurement strength 
is proportional to the drive on ai + , the drive strength 
therefore serves as an in-situ experimental knob for the 
phonon number variance. Selecting any values for the 
a,j + drive strengths, one may take a long time-average 
of both the output of a\ + and its squared value to re- 
spectively extract (n b ) c and (n b )1, where (u) c denotes 
an ensemble average of the operator u conditioned on 
the measurement record. Simultaneously, one may ob- 
tain (n^) c + A(n b ) c from a long time-average of the out- 
put of a 2+ , and use the information from ai+ to de- 
termine {n b ) ■ Thus one obtains sufficient information 
to determine (in a time-averaged sense) the steady-state 
conditional phonon number variance a\ = (ni) c — (n b ) 2 
for the selected values of the drive strengths. Repeat- 
ing this procedure for incrementally stronger values of 
the drive on ai+, one may observe the controlled col- 
lapse of the phonon number wavefunction. By relying 



on the ai+ measurement channel to collapse the wave- 
function, this protocol accommodates the fact that the 
experimentally observed g 2 is very weak [9]; sufficiently 
long time-averages of the output of a 2+ can yield the 
desired information regardless of the value of g 2 . 

Being able to collapse the wavefunction in this man- 
ner, however, implies certain parameter constraints. The 
study in [lj established two fundamental conditions for 
observing phonon number quantum jumps via continu- 
ous measurement of the output of ai+: both the n b mea- 
surement rate I\ and K\+ must be much greater than 
the mechanical Fock state decay rate. However, these 
conditions also apply for observation of phonon num- 
ber wavefunction collapse via continuous measurement 
because the jumps arise only when the phonon number 
variance is very small. The study in j4] considered the 
special case of a one-sided cavity where coupling to a 
thermal bath was the only source of mechanical dissipa- 
tion. In the more realistic case that we consider here of 
a two-sided cavity with continuous sideband cooling of 
the mechanical mode, the fundamental conditions may 
be expressed as r^ max \ k 1+ > j th [(n th + l)n b + n th (n b + 

1)] + (li,R + 12,r + lo P t)n b , where r^ max) denotes the 
maximum attainable value of T\, 7th is the mechanical 
dissipation rate to the mechanical mode thermal bath of 
average occupation n t h, Jopt is the mechanical dissipa- 
tion rate to the zero temperature optical cooling bath 
induced by the sideband cooling [11], and 7^ are due 
to the Raman scattering processes described above. Al- 
though 71. r is itself proportional to n\ + , which must be 
increased to increase Ti for each successive collapse incre- 
ment of the protocol, the total mechanical cooling rate 
7cooi = 7i, R + I2.R + lopt may be kept constant by si- 
multaneously reducing 7 opt via in-situ adjustment of the 
sideband cooling drive [IT] . 

The contraints on p( max ) anc i Kl + show that the fea- 
sibility of the scheme entails that n b be small. From de- 
tailed balance, n b = (j cooi n opt +7th«th)/(7cooi + 7th)- As 
the optical bath occupation n opt is very small, choosing 
7cooi ~ 7th^th yields n b w 1 and T^ max) ,K 1+ > 47 th n th . 
Thus, steady-state n b ss 1 can be achieved via continuous 
sideband cooling that is simultaneous with the collapse 
measurement and quantum jump measurement proto- 
cols without significantly increasing p( max ) beyond what 
would be required in the absence of continuous sideband 
cooling where the mechanical mode was instead passively 
cooled to n b ~ 1. Observing phonon number quantum 
jumps and wavefunction collapse with simultaneous side- 
band cooling may prove to be an experimentally more 
viable route than with passive cooling. 

The authors of [8] [10] derive the additional condition 
gi > for energy quantization detection by requir- 

ing that the phonon number measurement rate be greater 
than the mechanical Fock state decay rate due to the 
Raman process described above. However, because the 



4 



measurement plays the dual role of detecting the Fock 
state and also collapsing the distribution to create the 
Fock state, what is actually required is that the mea- 
surement rate be much greater than the Fock state de- 
cay rate. This was established in 0] and is reflected in 
the constraint on p( max ) above. The true requirement is 
therefore g\ 3> k + K-. 

Numerical demonstration. — In contrast to [8], we con- 
sider the case where the phonon measurement optical 
modes dj+ are strongly driven so that a J+ — > ctj + a'j + , 
where ctj are the steady-state background amplitudes 
of aj + and a'j + are the quantum fluctuations on top 
of Oj. We may therefore proceed in analogy to [4] 
to move to a displaced picture for modes a J+ and use 
the following stochastic master equation [T2] (SME) 
for the system density matrix p to treat the transmit- 
ted outputs of Oj+ as continuously observed via ho- 
modyne detection and the rest of the dissipative chan- 
nels as unobserved: dp = —^[H\,p]dt + (7th ("th + 
1) + jcoo\)'D[b}pdt + -f th n til V[b^}pdt + + 

K +M)'D[a J +\pdt + Y, 2 j =1 y /'n K j+,tdWj'H[aj + e~' l ^}p. Here 
we have dropped the prime from a'j + for convenience, 

Hi = -h*±{a\ + + a 1+ )n b - h^(a\ + + a 2 +)[nl + An b ) 
is the linearized interaction Hamiltonian in the displaced 
picture, \j = %Gj a ji ^[ c ]p = cpS — c i cp/2— pc^c/2 is the 
dissipation superoperator, H[c]p = cp+ pc) — Tr(co+ pc^) 
is the measurement superoperator, rj is the efficiency of 
the detectors, and dWj are independent Wiener incre- 
ments. Each optical mode a$+ has three dissipation 
channels at zero temperature with corresponding dissi- 
pation rates: reflected signal (Kj+. r ), transmitted sig- 
nal = K j+ ~ K j+,r)i and Raman scattering de- 
cay as discussed above. The mechanical mode 
b has four dissipation channels: thermal bath dissipa- 
tion (7th) and cooling dissipation (7 CO oi), consisting of 
sideband cooling and two Raman scattering channels. 
As explained above, 7 CO oi may be considered constant 
over the entire collapse measurement process. The pho- 
tocurrent homodyne measurement signals |12| are given 
by ij(t) = nKj+^aj+e'^ + a] + e^) c (t) + y/r)Kj +t i€j(t), 
where = dWj/dt. 

As per above, a% + must adiabatically follow the me- 
chanical mode. For computational simplicity we assume 
that 02+ does as well. From the quantum Langevin 
equations we find ai + = i p nh and a?+ = 

i p 2 (n? + Arih). Using this and the fact that the 

mechanical mode density matrix remains diagonal, we 
find 

dp n =7thra t h[np„-i - (n + l)p n ]dt 

+ (7th (nth + 1) + 7cooi) [(n + l)p n +i - n Pn] dt 

- 2^nT 2 {{n 2 + An) - (n 2 + An) c ) Pn dW 2 . (4) 
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FIG. 1. Collapse of phonon number distribution with increas- 
ing nb measurement strength (Ti) and constant n% +71& mea- 
surement strength (r2). For each value of Ti, the system is 
first allowed to reach steady-state, then (nt,) c , {nb)i, and {nl) c 
are obtained by computing (rib)c and (nf} c at each subse- 
quent time step and averaging for a sufficiently long time. The 
steady-state conditional variance decreases monotonically as 
a function of Fi, and becomes much less than the uncondi- 
tional variance when Ti 3> 47 t hMth, which is the regime in 
which quantum jumps arise in (n&) e (t) [!]■ The values of 
{nb) c indicate that the measurement processes shift the me- 
chanical mode out of perfect equilibrium with its effective 
thermal bath. Simulation parameters: T2 = 10 _5 (47 t hn t h), 
nth = 5 x 10 3 , 7cooi = 7thn t h, and A = l. 

2 

where T, = -, 3 , 3+,t , 2 , and p n is the occupation prob- 

ability of the nth mechanical Fock state. The photocur- 
rcnts become ii(t) = 2rjxi(nb) c (t) + y^/«r+~t^i (t) and 
i 2 (t) = 2r/X2("-b + An b)c{t) + v / T?«;2+ ) i6(i)- Here, Ti is 
the same phonon number measurement rate discussed in 
the previous section. Provided that A (defined above) 
is known, for fixed values of Tj sufficiently long time- 
averages of ii(t), and i 2 {t) yield the values of (n&) c , 
(rib) 2 , and (n|) c . We remind the reader that Tj are pro- 
portional to ctj , which can be adjusted in-situ by varying 
the optical drive strengths. 

Assuming T w 300 mK, n w 2tt x 1 MHz [<2 [13], 
and 7 th = 27r x 0.1 Hz [5], we find n th = 5 x 10 3 . As 
discussed above, we can set 7 CO oi = 7th«th so that nj, « 1. 
Arbitrarily setting A = 1, we assume n\ + ^> 47th«th so 
that Eq. Q is valid and we numerically integrate it for 
different values of Ti to produce the data shown in Fig. 

HI 

Conclusion and outlook. — This paper puts forth a new 
paradigm in which continuous weak measurement is used 
as an in-situ knob to control the quantum noise pro- 
file of a macroscopic degree of freedom, thus adding an 
additional capability to the toolbox of engineered quan- 
tum systems. The particular scheme presented is not far 
from experimental reality, as the latest iteration of the 
membrane-in-the-middle system in |13j shows impressive 
gains in the parameter space. A further order of mag- 
nitude increase in g\ and optimization of Ki± should 
achieve gi ^> 
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While in this paper the finite steady-state variance is 
the result of measurement-induced collapse being coun- 
teracted by a dissipative process, one can also imagine 
a scenario where the measurement-induced collapse in 
one variable may be counteracted by a continuous mea- 
surement on its conjugate variable. Also, in a many- 
body scenario changes in the quantum noise profile of 
degrees of freedom can lead to quantum phase transi- 
tions that dramatically change the bulk properties of the 
system, and it is interesting to consider the possibility 
of using the paradigm introduced in this paper to effect 
measurement-induced quantum phase transitions of en- 
gineered quantum metamaterials. 

A. A. G. acknowledges the support of the Australian 
Research Council Centre of Excellence for Engineered 
Quantum Systems and the UQRS/UQIRTA scholarships, 
and is grateful to Gerard J. Milburn for discussions and 
critiques of the manuscript. 



* [a.gangat@physics. uq .edu.au] 
[1] T. J. Kippenberg and K. J. Vahala, Science 321, 1172 
(2008); F. Marquardt and S. M. Girvin, Physics 2, 
40 (2009); M. Aspelmeyer et al, J. Opt. Soc. Am. B 
27, A189 (2010); G. J. Milburn and M. J. Woolley, 
Acta Physica Slovaca 61, 483 (2012); M. Aspelmeyer, 
P. Meystre, and K. Schwab, Phys. Today 65, 29 (2012); 



[2 
[3 

[4: 

[5 

[6 
[7 

[8 
[9 
[10 
[11 
[12 

[13 



M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, 
larXiv: 1303. 07331 

A. A. Clerk, F. Marquardt, and J. G. E. Harris, Phys. 
Rev. Lett. 104, 213603 (2010). 

A. M. Jayich, J. C. Sankey, B.M. Zwickl, C. Yang, J. D. 
Thompson, S. M. Girvin, A. A. Clerk, F. Marquardt, and 
J. G. E. Harris, New J. Phys. 10, 095008 (2008). 
A. A. Gangat, T. M. Stace, and G. J. Milburn, New J. 
Phys. 13, 043024 (2011). 

S. Mancini, V. Giovannetti, D. Vitali, and P. Tombesi, 
Phys. Rev. Lett. 88, 120401 (2002). 

Pikovski, I., M. Vanner, M. Aspelmeyer, M. Kim, and C. 
Brukner, Nature Phys. 8 (2012). 

J. D. Thompson, B. M. Zwickl, A. M. Jayich, F. Mar- 
quardt, S. M. Girvin, and J. G. E. Harris, Nature (Lon- 
don) 452, 72 (2008). 

M. Ludwig, A. H. Safavi-Naeini, O. Painter, and F. Mar- 
quardt, Phys. Rev. Lett. 109, 063601 (2012). 
J. C. Sankey, C. Yang, B. M. Zwickl, A. M. Jayich, and 
J. G. E. Harris, Nature Phys. 6, 707 (2010). 
H. Miao, S. Danilishin, T. Corbitt, and Y. Chen, Phys. 
Rev. Lett. 103, 100402 (2009). 

F. Marquardt, J. P. Chen, A. A. Clerk, and S. M. Girvin, 
Phys. Rev. Lett. 99, 093902 (2007). 

H. M. Wiseman and G. J. Milburn, Quantum Measure- 
ment and Control (Cambridge University Press, Cam- 
bridge, 2009). 

N. E. Flowers-Jacobs, S. W. Hoch, J. C. Sankey, A. 
Kashkanova, A. M. Jayich, C. Deutsch, J. Reichel, and 
J. G. E. Harris, Appl. Phys. Lett. 101, 221109 (2012). 



